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The hydrogen negative ion H“ is the simplest two-electron system that exists in nature. This 
system is not only important in astrophysics but it also serves as an ideal ground to study electron- 
electron correlations. The peculiar balance of the correlations between the two electrons with the 
interaction of electron-nucleus in H~ makes this system to have only two bound states, one being 
the ground state Is 2 1 S e and the other the doubly-excited metastable state 2p 2 3 P e embedded below 
the hydrogen n = 2 threshold. Here we report a calculation for the 2 p 2 3 P e state of H~ that 
yields the energy eigenvalue E = —0.125 355 451242 864 058 376 012 313 25(2), in atomic units. Our 
result substantially improves the best available result by 16 orders of magnitude. We further study 
the critical nuclear charge Z CI , the minimum value of nuclear charge Z that is required to bind a 
nucleus and two electrons. Our determination of Z CT for the 2 p 2 3 P e state of two-electron systems is 
= 0.994 781292 240 366 246 3(1), corresponding to 1/Z CI = 1.005 246 085 546 985 509 4(1), which 
improves the best published value of Z CI by about 10 orders of magnitude. We further investigate 
in a definitive way the unexplored regime of Z < Z CI using the method of complex scaling and 
establish precise shape resonance poles for the state of 2 p 2 i P e in the complex energy plane. 

PACS numbers: 31.15.ac 


Two-electron atomic systems that exist in nature in¬ 
clude the sequence H“, He, Li + ,..., which can be de¬ 
scribed by the following single Hamiltonian in nonrela- 
tivistic and infinite nuclear mass limit (in atomic units 
throughout) 


H = 



( 1 ) 


where iq and r 2 are the position vectors of electrons 1 
and 2 and ri 2 = |ri — r 2 1- There are two special bound 
states that always exist in the spectra of H for the above- 
mentioned sequence: one is the ground state Is 2 1 S e , and 
the other is 2 p 2 3 P e that is embedded between the n = 1 
and n = 2 one-electron ionization thresholds. These two 
states are truly bounded in the sense that they are sta¬ 
ble against autoionization. Moreover, for the first mem¬ 
ber in the two-electron sequence H _ , the state of Is 2 1 S e 
is the only bound state that exists below the n = 1 
ionization threshold, which has been proven mathemati¬ 
cally [Tj; and the state of 2 p 2 3 P e is the only bound state 
embedded in the continuum that has been firmly estab¬ 
lished numerically so far (2H81- However, as Z further 
reduces from Z = 1, these two states will eventually be¬ 
come unstable. For a given state, the minimum value of 
Z that the Hamiltonian can still keep this state to be 
bounded is called the critical nuclear charge Z CI of the 
state of interest. A precise determination of the critical 
nuclear charge is essential in understanding the analyt- 
icity of the function E[Z) or E{ A) in the complex Z or 
A plane, where A = 1/Z. In fact, for the ground state, 
A s = 1 / Z CT is the radius of convergence for E{ A) about 


A = 0 in the complex A plane, where E{ A) is analytic 
inside the circle |A| < A s but has an essential singularity 
at A = A s ~ 1.097660 833 [M2- 

For H _ , the energy eigenvalue of Is 2 3 S e has been cal¬ 
culated to high precision m M- The corresponding 
critical nuclear charge has been investigated by several 
authors, with the most precise calculation of Estienne et 
al. [12 that sets a firm benchmark of 16-digit accuracy. 
Atomic properties near and at the critical nuclear charge 
of Is 2 1 S e have recently been studied by Grabowski and 
Burke m- When Z < Z CI , the situation becomes more 
complicate and subtle, as the bound state now crosses 
over the nearest ionization threshold and turns into a 
shape resonance. Since a shape resonance lies in a scat¬ 
tering continuum, its energy, usually complex, cannot be 
defined directly from the original Hamiltonian |l]), but 
can be determined from its proper analytical continua¬ 
tion Em. These resonances have been explored by Dubau 
and Ivanov m, and Sergeev and Kais |T6] using the 
method of complex scaling. Comparing to the ground 
state, the analytic properties and resonant structure of 
the 2 p 2 3 P e state is much poorly understood. There are 
only two calculations available in the literature, to the 
best of our knowledge: one by Brandas and Goscinski El 
who give l/^cr = 1.0048 and one by Sergeev and Kais [T8] 
who give 1 / Z cr = 1.00524608. Also to our knowledge, 
no studies have been reported to date concerning res¬ 
onances in the region Z < Z CT for the pure Coulomb 
case, notwithstanding that an investigation of shape res¬ 
onances for the screened Coulomb case has been reported 
in Ref. [15] - 
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The purpose of the present Letter is threefold. First, 
we present a substantially improved calculation on the 
energy eigenvalue of the 2 p 2 3 P e state in H~. Second, we 
establish a much more definitive determination for the 
critical nuclear charge of this state. Finally we explore, 
for the first time, the unexplored region where Z < Z CT 
and report the 2 p 2 3 P e resonant energies and widths using 
the method of complex scaling. In all these calculations, 
we employ Hylleraas bases to describe the correlations 
between two electrons. 

We use the following basis set in Hylleraas coordinates 

r i r 2 r i 2 exp(-ar! - /3r 2 ) (f i, r 2 ) , (2) 

where a > t\ and b > t 2 , 

h,h)= Y, (hm. 1 £ 2 m. 2 \LM)Y^(h)Y^(r 2 m 

m\m2 

is the common eigenfunction of the total angular mo¬ 
mentum squared L 2 , the z component L z , and the parity 
operator with the corresponding eigenvalues L(L + 1), 
M, and (—l) £l+ ^ 2 , respectively. The size of basis set is 
determined by including all terms with a + b + c < f 1 
with S2 being an integer. For the 2 p 2 3 P e state of H _ , we 
choose i\ = l 2 = L = 1. We include five blocks in the 
basis set each with different (a,/?). In order to maintain 
maximum numerical stability of the basis set, we only in¬ 
clude terms with a < b. Starting from the second block, 
all terms with a = b are also excluded. Table U shows 
a convergence study of the 2 p 2 3 P e energy eigenvalue as 
the size of basis set increases. In the table, i?(H) is the 
ratio of two successive differences that can be used as a 
measure for the rate of convergence. We thus stop the 
calculation where i?.(fi) ~ 1. The achieved accuracy of 
the extrapolated value of the energy is about 2 parts in 
10 28 . Comparison with available published results indi¬ 
cates that the present calculation has dramatically im¬ 
proved the previous best result of Kar and Ho 0 by 
about 16 orders of magnitude. 

We now turn to the problem of finding the criti¬ 
cal nuclear charge for the 2 p 2 3 P e state. The general 
condition for this state being bound below the hydro- 
genic n = 2 threshold is that E(Z) < —Z 2 /(2n 2 ), 
or e(Z) = -E(Z)/Z 2 > 1/(2 n 2 ) = 0.125. The ba¬ 
sis set used in the variational calculation of the criti¬ 
cal nuclear charge is similar to that of the ground state 
except in the case of having five blocks in the basis, 
we only use two blocks. The nonlinear parameters for 
the first block is approximately (a,/?) ~ (0.5,0.01), 
and for the second block (a,/?) ~ (0.7, 0.5). Table |ll| 
lists the scaled energy e(Z) for a series of calculations 
of different values of Z. We use the arbitrary preci¬ 
sion software package MPFUN, developed by Bailey [2TI| . 
The precision we set is 80 decimal places. Our re¬ 
sults in Table [II] clearly shows that the critical nuclear 
charge Z CI must be between 0.994 781292 240 366 246 0 
and 0.994 781 292 240 366 246 5. We thus give an extrapo¬ 
lated value of Z CI to be 0.994 781292 240 366 246 3(1) that 


contains at least 18 significant digits after the decimal 
point. Our result is not only in agreement with the value 
of Sergeev and Kais [H] but also more precise than theirs 
by 10 more significant figures. 

In order to explore the region of Z < Z CI , we ap¬ 
ply the method of complex scaling [2T], by which the 
inter-particle coordinates are transformed according to 
r —>• rexp(*0) with the phase angle 9 being real and 
positive. The Hamiltonian |l]) is then transformed as 
H —> H(9) = T exp(—2 id) + V exp(— i&), where T and 
V are the kinetic and potential operators respectively. 
A diagonalization of the transformed Hamiltonian yields 
complex energy eigenvalues of the form E res = E r — iT /2, 
where E r is the resonance position and T the resonance 
width. In search of a resonance, a resonance pole is the 
one that exhibits the most stabilized behavior of the en¬ 
ergy with respect to changes of 9 and nonlinear parame¬ 
ters in the basis set. Table Hill shows numerical results for 
some selected complex resonance energies with different 
Z < Z CT . It should be mentioned that we have used two 
different types of basis sets: the Hylleraas and the con¬ 
figuration interaction (Cl) and the results are consistent 
up to 4 ~ 5 digits for the shape resonances in Table EH 

Figure [l] is the plot of resonance energy versus 1/Z 
for this shape resonance. It is seen that the shape reso¬ 
nance starts to appear just above the two-body H (n = 2) 
threshold as Z is first decreased ( 1/Z is increased) be¬ 
low Z CT . As Z is decreased further, the shape resonance 
moves away from the two-body H (n = 2) threshold. 
When 1/Z is increased further to near the region around 
1.8, the shape resonance starts to approach the H (n = 3) 
threshold and it even crosses over the H (n = 3) thresh¬ 
old when 1/Z is increased to around 1.9. We stop at 
1/Z ~ 2.0, because as the resonance is approaching the 
total ionization limit E = 0, it is getting more difficult 
to obtain convergence in resonance calculations. 

Figure [2] is the plot of the half-width T/ 2 versus 1/Z 
in logarithmic scale. It is seen that the width increases 
very rapidly when the bound state first turns to a shape 
resonance as its energy crosses over the H (n = 2) thresh¬ 
old. The parent state of this shape resonance is hence the 
excited 2 p state of the two-body atom, together with the 
other electron having a p- wave character (two spin-half 
p-wave electrons are coupled to form a 3 P e state). Now 
when the resonance crosses the n = 2 threshold from 
a bound state to a shape resonance at Z cr , the p- wave 
electron will experience a repulsive angular momentum 
barrier, through which the electron can tunnel out, re¬ 
sulting in a cross-over from a bound state to a shape 
resonance. Furthermore, as the shape resonance moves 
away from the threshold when Z is decreased further, the 
thickness of the potential barrier, through which the elec¬ 
tron tunnels out, becomes narrower, leading to a shorter 
lifetime for the autoionization process, hence broadening 
the width, a consequence of the uncertainty principle. 
The physics of such a phenomenon is similar to the elec¬ 
tron tunneling in field ionization when an atom is placed 
under an external electric field. 
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One can see from Fig. [2] that F/2 reaches a maxi¬ 
mum value near 1 jZ ~ 1.7, and then it starts to de¬ 
crease somewhat as \/Z further increases. This can be 
explained as follows. As shown in Fig. [TJ when Z is 
decreased below Z CI , a bound state becomes a shape res¬ 
onance lying above the H (n = 2) threshold. As Z is 
decreased further, the shape resonance is moving away 
from the H (n = 2) threshold and approaching the H 
(n = 3) threshold. As such, the resonance is now gain¬ 
ing some Feshbach-type character. In other words, at 
this range of 1 /Z, the resonance is a mixture of a shape 
resonance, with the H (n = 2) as the parent state, and 
a Feshbach resonance, with the H (n = 3) as the par¬ 
ent state. The autoionization process is a combination 
of tunneling effect, as described earlier, and Feshbach- 
type configuration-interaction effect as one electron de¬ 
caying from the n = 3 shell to the n = 2 shell, with 
the other electron being autoionized away. Usually, for 
a Feshbach resonance the width would decrease as it is 
moving toward its parent (n = 3, in this case) threshold. 
So for our special case, the total width is a competition 


between the “increasing character” due to tunneling ef¬ 
fect and the “decreasing character” due to the Feshbach 
configuration-interaction effect. Apparently, as the reso¬ 
nance approaches the H (n = 3) threshold, the Feshbach 
resonance effect would become more pronounced, result¬ 
ing in the overall decreasing trend for the total width 
around this 1 jZ region. 

In conclusion, we have established two new benchmark 
results for the 2p 2 3 P e state, i. e, the energy eigenvalue of 
H _ and the critical nuclear charge for the two-electron 
isoelectronic sequence. We have also explored in a defini¬ 
tive way the resonance poles as the nuclear charge is 
below its critical value. Our results will be valuable in 
studying the analytical structure of E{ A) in the complex 
A plane. 
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TABLE I: Convergence of the nonrelativistic energy for the 
2 p 2 i P e state of H“ with infinite nuclear mass. In atomic 
units. 


n 

No. of terms 

E{Q) 

R{fl) 

14 

1064 

-0.125 355 451 242 864 057 032 738 135 921 31 


15 

1316 

-0.125 355 451 242 864 058 265 245 521423 53 


16 

1604 

-0.125 355 451 242 864 058 368 324 462 549 58 

11.95 

17 

1932 

-0.125 355 451 242 864 058 375 223 510 548 96 

14.94 

18 

2301 

-0.125 355 451 242 864 058 375 939 082 805 06 

9.64 

19 

2715 

-0.125 355 451 242 864 058 376 004 981 517 50 

10.85 

20 

3175 

-0.125 355 451 242 864 058 376 011 372 600 17 

10.31 

21 

3685 

-0.125 355 451 242 864 058 376 012 166 845 94 

8.04 

22 

4246 

-0.125 355 451 242 864 058 376 012 296 758 16 

6.11 

23 

4862 

-0.125 355 451 242 864 058 376 012 311 685 94 

8.70 

24 

5534 

-0.125 355 451 242 864 058 376 012 312 844 98 

12.87 

25 

6266 

-0.125 355 451 242 864 058 376 012 313 076 82 

4.99 

26 

7059 

-0.125 355 451 242 864 058 376 012 313 169 71 

2.49 

27 

7917 

-0.125 355 451 242 864 058 376 012 313 202 85 

2.80 

28 

8841 

-0.125 355 451 242 864 058 376 012 313 222 51 

1.68 

29 

9835 

-0.125 355 451 242 864 058 376 012 313 231 80 

2.11 

30 

10900 

-0.125 355 451 242 864 058 376 012 313 237 21 

1.71 

Extrap. 


-0.125 355 451 242 864 058 376 012 313 25(2) 


Aashamar [2] (1970) 


-0.125 325 


Drake [3] (1970) 


-0.125 350 


Bhatia |4j (1970) 


-0.125 354 7051 


Jauregui and Bunge 0 (1979) 


-0.125 354 716 6 




-0.125 355 08(10) (extrap.) 


Banyard et al. [Bj (1992) 


-0.125 353 6 


Bylicki and Bednarz [7j (2003) 


-0.125 355 45124 




-0.125 355 453 06 (extrap.) 


Kar and Ho [H] (2009) 


-0.125 355 451242 








TABLE II: Scaled energy eigenvalue e(Z) = —E(Z)/Z 2 for 
the state of 2 p 2 2 P e as a function of Z. The size of basis set 
in each calculation is 11900. In atomic units. 
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Z 

0.994 781292 240 370 000 0 
0.994 781292 240 367 000 0 
0.994 781292 240 366 400 0 
0.994 781292 240 366 250 0 
0.994 781292 240 366 247 0 
0.994 781292 240 366 246 5 
0.994 781292 240 366 246 0 
Extrapolation: 

0.994 781 292 240 366 246 3(1) 
0.994 78129° 


_ <Z) _ 

125 000 000 000 000 228 829 28 
125 000 000 000 000 048 404 64 
125 000 000 000 000 009 865 61 
125 000 000 000 000 000 230 86 
125 000 000 000 000 000 038 16 
125 000 000 000 000 000 006 05 
124 999 999 999 999 999 973 93 


0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 


“Ref. EE]. 


TABLE III: Selected shape resonance poles for the state 
2 p 2 i P e for various Z values when Z < Z cr . In the table, 
a[6] = a x 10 6 . In atomic units. 

Z W r F 72 

0.993 —1.231466[—1] 4.9685[-6] 

0.990 —1.222194[—1] 2.8275[—5] 

0.985 —1.206725[—1] 8.7928[-5] 

0.980 —1.191255[— 1] 1.6185[—4] 

0.970 —1.160375[— 1] 3.3449[-4] 

0.960 —1.129645[— 1] 5.2682[-4] 

0.940 —1.068834[—1] 9.4119[-4] 

0.920 —1.009096[—1] 1.3716[-3] 

0.900 —9.505911 [—2] 1.8031[—3] 

0.850 —8.103785[—2] 2.8364[-3] 

0.800 —6.795992[—2] 3.7552[-3] 

0.700 —4.486793[—2] 5.1245[—3] 

0.600 —2.607694[—2] 5.7485[-3] 

0.500 —1.170528[—2] 5.5065[-3] 



















